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ABSTRACT 



We present the results of a spectroscopic study of the Fe Ka emission of the persistent neutron-star atoll low-mass X-ray binary 
and type I X-ray burster GX 3 + 1 with the EPIC-PN on board XMM-Newton. The source shows a flux modulation over several years 
and we observed it during its fainter phase, which corresponds to an X-ray luminosity of L x ~ 10 37 ergs s~ l . When fitted with a 
two-component model, the X-ray spectrum shows broad residuals at ~ 6-7 keV that can be ascribed to an iron K a fluorescence 
line. In addition, lower energy features are observed at ~ 3.3 keV, ~ 3.9 keV and might originate from Ar XVIII and Ca XIX. The 
broad iron line feature is well fitted with a relativistically smeared profile. This result is robust against possible systematics caused 
by instrumental pile-up effects. Assuming that the line is produced by reflection from the inner accretion disk, we infer an inner disk 
radius of ~ 25R g and a disk inclination of 35° < i < 44°. 
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1. Introduction 

Spectra of accreting compact objects, such as active galactic 
nuclei and X-ray binary systems containing stellar mass black 
holes or weakly magnetized neutron stars (NS), often exhibit 
broad emission lines (FWHM up to ~ 1 keV) at E ~ 6.4 - 
6.97 keV, identified with K a radiative transitions of iron at dif- 
ferent ionization states. These broad lines are thought to orig- 
inate from reprocessed emission that is relativistically smeared, 
in the inner parts of the accretion disk illuminated by the primary 
Comptonized spectrum. This process produces different spectral 
signatures, including a reflection hump that peaks between 20- 
30 ke V, emission lines, and absorption edges (see iFabian et al.1 
120001 for a review). Owing to their relatively high abundance and 
fluorescence yield, the iron lines are the best suited features to di- 
agnose the accretion flows close to the compact object, because 
the predicted line profile depen ds on the geometry - inclination 
and inner radius - of the system (IFabian et al. 1989). For NS sys- 
tems, the inner disk radius sets an upper limit to the radius of the 
star, which can constrain the NS equation of state. 

iPiraino et al. I d2000h first suggested a relativistic origin of the 
broad iron line observed in the NS system 4U 1728-34. Recently, 
thanks to the improved spectral capability of XMM-Newton and 
Suzaku asymmetric and smeared broad Fe K lines have been ob- 
serve d in about 10 NS l o w-mass X-ray binaries (|Di Salvo et aD 
2009t iD'Ai eTail 120091: llaria et al l l2009t iPande et al l 120081 



Bhattacharvva & Strohmaverl 120071 ICackett et al.l 120081 Boiol 
and were interpreted as relativistic reflection features from the 
inner part o f the accretion dis k. As shown in the self-consistent 
approach of lD'Ai et al.l (1201 Oh . the broad features other than that 
of iron in the spectrum also support the "relativistic origin" sce- 



nario. However, the robustness of t hese resul t s with respect to 
pile-up effects has been questioned. iNg et al.l d2010l) . based on 
archival XMM observations, studied the Fe K emission of 16 
NS systems, and concluded that the line profiles show no ev- 
idence of asymmetry w hen spectra a r e extr acted from pile-up 
free regions. Conversely. iMiller et al.l (1201 Ol) conducted simula- 
tions for CCD based detectors to assess the impact of photon 
pile-up on relativistic disk line and continuum spectra, conclud- 
ing that relativistic disk spectroscopy is generally robust against 
pile-up when the effect is modest. 

Although the "relativistic origin" is a widely accepted inter- 
pretation, alternative physical mechanisms to explain the line 
broadening are still debated. A broad iron line can originate 
from thermal Comptonization of line photons in th e so-called 
accretion disk corona (ADC) dBrandt & Matt! 1 1994 . or can be 
caused by non-thermal Com ptonization in strong outflow winds 
dLaurent & Titarchu"kll2007h . The predicted line profiles differ 
significantly from the relativistic disk line profiles. 

In this Letter, we report detection of a broad and asymmetric 
iron line in a spectrum of the neutron star LMXB GX 3 + 1 with 
high statistics and excellent energy resolution obtained with the 
EPIC-pn camera on board XMM-Newton. We discuss in detail 
instrumental pile-up effects and show that the relativistic profile 
of the line is robust against biases from photon pile-up. 



2. Observation and data analysis 

GX3+1 is a typical bright (-F2-10 keV ~ 150 - 400 mCrab) atoll 
source with peculiar bursting behavior. It exhibits type -I bursts 
dMakishima et alJll983b iKuulkers & van der Klisll2000h and su- 
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perbursts, thought to arise from car bon shell flashes in the layers 
below the surface (lKuuikersll2002l) . The spectrum of the source 
observed wi th the BeppoSAX narrow field instruments has been 
modeled by lOosterbroek et al.l (1200 ll) with a disk-blackbody, a 
Comptonized component, and absorption by interstellar medium. 
Interestingly, a broad iron line was clearly detected but only 
marginally discussed by the authors. 

GX3+1 was observed on September 1, 2010 by XMM- 
Newton and simultaneously with RXTE for about 56 ks. 

The XMM-Newton Observatory dJansen et al.l 1200 ll) in- 
cludes three 1500 cm 2 X-ray telescopes each with a European 
Photon Imaging Camera (EPIC, 0.1-15 keV) at th e focus. Two 
of the EPIC imaging spectro meters use MPS CCD s dTurner et al.l 
1200 lb and one uses pn CCDs dStriider et al.l2001l) . Data products 
were reduced using the science analysis software (SAS) version 
1 1 .0. Owing to the expected high count rate from the source, the 
telemetry of EPIC MOS cameras was allocated to the EPIC pn 
camera, therefore no data from the MOS cameras are available. 
We did not analyze the RGS data since they do not cover the 
Fe K energy band. 

The EPIC pn was operated in timing mode for the entire ob- 
servation. In this mode only one CCD chip is operated and the 
data are collapsed into a one-dimensional row (4-'A) and read out 
at high speed, the second dimension being replaced by timing 
information. This allows a time resolution of 30 /-is. We used 
the SAS task epfast on the event files to correct for the ef- 
fect of charge transfer inefficiency (CTI), which has been seen 
in the EPIC pn timing mode at high count ratefl Ancillary re- 
sponse files were generated using the SAS task arfgen follow- 
ing the recommendations in the XMM-Newton SAS User guided 
for piled-up observations in timing mode whenever applicable. 
Response matrices were generated using the SAS task rmfgen. 

Since the average EPIC pn count rate was close to the 
600 counts s level at which pile-up effects start to be impor- 
tant, we investigated in detail if there was any pile-up before 
extracting the spectra. We used the SAS task epatplot, which 
uses the relative ratios of single- and double-pixel events as a 
diagnostic tool. These ratios deviate from standard values when 
there is significant pile-up in the pn camera timing mode. 

The source coordinates fell in the center of column 37 of 
the CCD. Source events were first extracted from a 53.3" (13 
colums) wide box centered on the source position; we call this 
Region 1 . We also defined regions excluding the central columns. 
Region 2 has the central column (column 37) removed. Region 3 
has the three central columns removed (columns 36, 37, and 38). 
The relative ratios of single- and double-pixel events exhibits an 
energy-dependent deviation between the data and the model at 
E > 4 keV when Region 1, including all source events, is used. 
The effects of pile-up are reduced for Region 2, with the central 
column removed. After removal of the three central columns, 
Region 3, any deviation in the ratios disappears, showing that 
pile-up has no significant effect. Spectra were extracted for each 
of the three regions, selecting only events with PATTERN < 4 
(singles and doubles) and FLAG=0. The energy channels were 
grouped by a factor of four to avoid oversampling the instrumen- 
tal energy resolution. 

3. Results 

To model the spectrum, we used the spectral analysis package 
XSPEC v. 12.7.0. To avoid the systematic effects discussed in 

1 http://xmm.esac.esa.int/external/xmm_calibration 

2 http://xmm.esac.esa.int 



iD'Ai et aT] (1201 Ol) that were recently addressed in the EPIC-pn 
calibration notqj we removed the energy bins below 2.4 keV 
and those above 9 keV. We then fitted t he continuum with 
a model widely us e d for atoll LMXBs dPiraino et al ] I2007t 
iDi Salvo et all 120091; iBarret & Olivell2002l) : a photoelectric ab- 
sorbed two-component model (Model 1) containing a blackbody 
and a th ermal Comptoni zed component described by the model 
comptt dTitarchukl 1 1 9941) . For the photoelectric absorption, we 
used the XSPEC model phabs. 

This model gave an unacceptable fit (^ 2 /DoF = 2621/322) 
because there were several unmodeled features across the entire 
energy range, the most prominent being a broad profile with 
residuals from 5.5 to 7.3 keV and a maximum deviation of 
~ 14<x from the best-fit co ntinuum (see F i g. [1] b lack points). 
Following the approach of |Di Salvo et al . (2009), we added 
to the continuum three Gaussian emission lines at ~3.3 keV 
(Ar XVIII), -3.9 keV (Ca XIX), and ~ 6.6 keV (Fe XXIV 
and/or Fe XXV), respectively. This greatly improved the fit to 
;^ 2 /DoF = 557/313. All these lines are significantly broader 
than the pn energy resolution, with Gaussian a of 140 eV for 
Ar, 180 eV for Ca, and 290 eV for Fe. The fit was additionally 
improved by adding an absorption edge at 8.40 keV resulting in 
;r/DoF = 521/311. 

Because the residuals of the Fe line profile, shown in Fig. [TJ 
suggest an asymmetric shape, w e replaced the Gau ssian at 
6.6 keV with a diskline profile dFabian et al.l 1 19891) . obtain- 
ing again a significant improvement of the fit with ^- 2 /DoF = 
45 1/311. Adding two parameters improved the fit by A^ 2 = 100, 
corresponding to an F-test probability of chance improvement 
of 3 x 10 ). However, the comptt parameters are poorly con- 
strained and the high value of the optical depth, near 10, sug- 
gests that the Comptonized component is close to a Wien spec- 
trum. Therefore, we used a second two-component continuum 
model containing a diskbb component and a simple blackbody 
component (Model 2). For this model all continuum parameters 
are well constrained. 

As for Model 1, using a diskline profile for the iron emission 
line added to Model 2 significantly improves the^- 2 /DoF with 
respect to a Gaussian from 549/315 to 441/313. The A^ 2 = 
108 for the addition of two parameters corresponds to an F-test 
probability of chance improvement of 1.3 x 10 . The diskline 
parameters are well constrained, except for the outer disk radius, 
which we fixed after some trials at R o „,=4500 (GM/c 2 ). The best- 
fit parameters are reported in Table QJand the unfolded spectrum 
and residuals are shown in the upper panel of Fig. |2] 

The best-fit value obtained for the energy of the iron line, 
Eune, suggests an ionization parameter of log £ 2.5 -2.7, with £ 
defined as Lx/{nr 2 ), where Lx is the ionizing X-ray luminosity, n 
the density in the reflector, and r 2 is the distance from the central 
source of emission. 

In Column 2 of Tab. [TJ we report the best-fit parameters of the 
spectrum of GX 3+1 for data extracted in Region 2 (without the 
central column 37). All spectral parameters are compatible with 
those obtained for data extracted from Region 1. The diskline 
gives again a better fit, with a A^ 2 = 73 for the addition of two 
parameters, corresponding to an F-test probability of chance im- 
provement of 2 x 10 -12 . 

In the bottom panel of the Fig. [2 we show the unfolded pn 
spectrum of GX3+1 with individual components and residuals 
in units of cr for data extracted from Region 3 (without the cen- 
tral columns 36, 37, 38). The best-fit parameters are reported in 
column 3 of Tab. [TJ Leaving all parameters free and using the 

3 http://xmm2.esac.esa.int/docs/documents/CAL-TN-0083.pdf 
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Energy (keV) 



Fig. 1. Residuals in units of cr in the pn energy range used with 
respect to the best-fitting continuum model, composed of an ab- 
sorbed (phabs) diskbb and a blackbody component for the spec- 
trum extracted in Region 1 (black), Region 2 (red) and Region 3 
(blue). 



2.4-9 keV energy range, the Nh value and the diskbb normal- 
ization slightly differ from those obtained for Regions 1 and 2, 
while the line parameters remain consistent. However, if we use 
data up to 11 keV, the Nh value and the diskbb normalization 
agree with those obtained for the other regions, since this spec- 
trum is not affected by pile-up even at high energies, . Therefore, 
we fixed Nh to the value obtained with data up to 11 keV, and 
again fitted the spectrum in the energy range 2.4-9 keV. All pa- 
rameters are now consistent with those obtained for Region 1 , as 
shown in Tab. Q] column 3. The diskline model improves the fit, 
with a Ax 2 = 36 for the addition of two parameters and an F-test 
chance improvement probability of 2 x 1CT 7 . 

As shown in Fig.Q] where we show the residuals in units of 
cr with respect to the best-fitting continuum model, for the spec- 
trum extracted in Region 1 (black), Region 2 (red) and Region 3 
(blue), the line shape does not strongly depend on the photon 
pile-up fraction. To check if the low-energy lines are affected 
by Doppler and relativistic distortion, we replaced the Gaussian 
lines in Model 2 with disk lines, with all the parameters fixed at 
those of the iron diskline except the line centroid energy and the 
normalization (Model 3). We obtained an equally good fit with 
the elimination of two parameters, chi 2 /DoF = 369/316 (see the 
rightmost column of Table[T). 



4. Discussion 

We reported the detection of broad emission lines at ~3.3 keV, 
-3.9 keV and -6.6 keV, corresponding to Ar XVIII, Ca XIX, 
and most likely a blend of Fe XXIV and Fe XXV, in high- 
quality XMM-Newton X-ray spectra of the LMXB and atoll 
source GX 3+1 obtained when the source was in a relatively low- 
intensity state. A b road, cr ~ 1.0 keV, i ron li ne feature was pre- 
viously reported by lOosterbroek et al. (2001), who simply men- 
tioned that the line might originate from reflection from an ac- 
cretion disk or from Compton broadening in the external parts 
of a ~3 keV corona. We showed that the disk-reflection scenario 
is consistent with the high-quality data from the XMM-Newton 
EPIC-PN. The iron line's broadening and shape are therefore af- 
fected by Doppler and relativistic effects at the inner edge of the 
accretion disk. In addit ion, Compton broaden ing is expected due 
to the hot ionized disk (Ross & Fabianl200"l . Such a component 
does not affect the skewness of the line, however, and indeed we 
have verified that even including Compton broadening, relativis- 
tic smearing is required with high statistical significance. 



a. 
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Fig. 2. Top: unfolded pn spectrum and individual model com- 
ponents for Model 2 using data from Region 1 which includes 
piled-up events. Bottom: the same for Model 2, with Nh fixed 
as described in the text, for Region 3, which excludes piled-up 
events. Residuals in units of cr are shown in each panel. 



In our study of the spectral properties, we carefully inves- 
tigated effects caused by instrumental photon pile-up. We con- 
clude that the line shape does not significantly depend on the 
photon pile-up fraction. Our results agree with the conclusion of 
Miller et al. (2010) that relativistic disk spectroscopy is generally 
robust against pile-up when this effect is modest. We measured 
an iron line rest-frame energy, E ~ 6.61 keV that is compatible 
with line emission from Fe XXIV and Fe XXV, showing that 
the reflecting matter is ionized and leading to an estimate for the 
ionization parameter of log £ 2.4 - 2.7. The line appears to be 
produced at radii extending inward to Rj„ = (25 + 4-)R g , where 
R g = GM/c 2 is the gravitational radius. Because R g ~ 1.5 km for 
a 1M neutron star, the inner radius is then /?,„ = (50+8)oti.4 km, 
where mi a is the neutron star mass in units of 1.4 M . For rea- 
sonable values of the neutron star mass, our results indicates that 
the disk is truncated quite far from the NS surface. The apparent 
inner radius of the soft disk blackbody emission gives a radius 
of 8 + 2 km, calculated for a d istance to the source of 4.5 kpc 
dKuulkers & van der Kli si [20001) . Note that this value is thought 
to underestimate the tr ue inner radius by a factor of more than 
2 dMerloni et al.ll200"ol) . The power law dependence of the disk 
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Table 1. The best-fit parameters for spectra of GX 3 + 1 observed 
with the XMM-Newton EPIC-pn using different models and data 
extraction. The continuum used in both models reported here, 
Models 2 and 3, consists of the sum of a blackbody and a mul- 
ticolor disk blackbody. For Model 2 we report results obtained 
using data extracted from three different regions. In Model 3, 
the low-energy lines are modeled with diskline profiles, with all 
parameters fixed at those of the iron diskline except the line cen- 
troid energy and the normalization. The blackbody luminosity is 
given in units of L-yijD^, where L37 is the bolometric luminos- 
ity in units of 10 37 ergs/s and D10 the distance to the source in 
units of 10 kpc. The diskbb normalization is given in units of 
(R-,,,1 D\q) 2 cos6 where /?,„ is an apparent disk radius in km and 
6 the angle of the disk. The line intensity 7/,„ e is given in units 
ph crrT 2 s . is referred to the switch from a Gaussian to 
a diskline to model the iron line, and the Prob the relative 
probability of chance improvement given by the F-test. Flux is 
calculated in the 2-10 keV band, and given in units erg cirT 2 
s -1 . Uncertainties are given at a 90% confidence level. 



Parameter 




Model 2 




Model 3 




Region 1 


Region 2 


Region 3 


Region 3 


/V H [10 22 cnr 2 ] 


1.7 ±0.1 


1.7 ±0.1 


I. & fixed 


I. & fixed 


*r bb [keV] 


1.70 ±0.02 


1.68 ±0.02 


1.62 ±0.01 


1.60 ±0.01 


^bb 


4.4 ±0.1 


4.4 ±0.1 


4.4 ±0.1 


4.45 ± 0.05 


ir dbb [keV] 


0.93 ± 0.03 


0.92 ± 0.04 


0.89 ± 0.02 


0.85 ± 0.02 




228 ± 37 


234+50 

-41 


278+g 


32 1 +24 

J -25 


E lim (keV) 


3.37 ± 0.02 


3.38 ± 0.02 


3.38 ± 0.03 


3.32 ± 0.03 


cr[keV] 


0.12 ±0.03 


0.11 ±0.03 


0.09 ± 0.04 




htm (10- 4 ) 


9±2 


8±2 


8±2 


13 ±3 


EqW (eV) 


5± 1 


5±1 


5 ± 1 


8±2 


E,„„ (keV) 


3.95 ± 0.02 


3.95 ± 0.02 


3.94 ± 0.02 


3.91 ±0.03 


cr[keV] 


0.08 ± 0.03 


0.06 ± 0.03 


0.07 ± 0.04 




hne (l<r 4 ) 


6± 1 


5±1 


6±2 


12 ±2 


EqW (eV) 


5 ±3 


4±3 


5±3 


10 ±2 


E liM (keV) 


6.62 ± 0.02 


6.60 ± 0.02 


6.60 ± 0.03 


6.61 ±0.02 


hne (10- 3 ) 


2.4 ± 0.2 


2.6 ±0.2 


9 q+0.5 


2.6 ±0.2 


EW[eV] 


57 ±5 


63 ±5 


70 ±6 


64 ±5 


Betor 


-2.61 ±0.03 


-2.62 ± 0.07 


-2.6 ±0.1 


-2.54 ± 0.09 


R m (CM/c 2 ) 


25 ±4 


25+5 


28 -?2 


17+7 
^'-14 


Incl (deg) 


37 ±2 


39 ±2 




39 +4 


4r 


108 


73 


36 






i.3 x nr 15 


2x ur 12 


2 x 10~ 7 




Flux (10-' ) 


4.53 ± 0.05 


4.49 ± 0.05 


4.48 ± 0.06 


4.48 ± 0.06 


X 1 (dof) 


441 (313) 


389 (313) 


352 (314) 


369 (316) 



to present more extensive spectral-timing studies based on the 
simultaneous RXTE observation and on the XMM RGS in the 
near future. 
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emissivity, r , is parametrized by the index betor, which is 
-(2.61 ± 0.03), indicating that the di sk emission is dom inated by 
irradiation by a central source (see iFabian et al.lll989l) . Finally, 
the disk inclination with respect to the line of sight is also well 
determined, i = (37 ± 2) degrees. This agrees well with the inter- 
mediate inclinations, between 20 and 65 degrees, found in other 
sources, as indeed expected from selection criteria. 

Signatures of a disk reflection are also the discrete features 
present at soft X-rays from Ar XVIII, and Ca XIX, respectively. 
We note that the scenario emerging from this work on GX 3 + 1 is 
coherent with results obtained for other sou rces and is in particu - 
lar very similar t o the case of 4U 1705-44 dDi Salvo et al.ll2009b 
iD'Ai et al.l[2QlQh . 

We focused on the robust detection of the relativistic iron line 
in the source, taking into account systematic instrumental effects 
caused by pile-up. Additional studies based on more physical 
self-consistent reflection models are in progress. We also aim 
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